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Abstract. In this paper we apply the reverberation method to determine the mass of the black hole in 3C 273 from the Lya and 
Civ emission lines using archival IUE observations. Following the standard assumptions of the method, we find a maximum- 
likelihood estimate of the mass of 6.59 10 9 M Q , with a lcr confidence interval 5.69 - 8.27 10 9 M Q . This estimate is more than 
one order of magnitude larger than that obtained in a previous study using Balmer lines. We reanalyze the optical data and show 
that the method applied to the Ha, H/3, and Hy Balmer lines produce mass estimates lower by a factor 2.5, but already much 
larger than the previous estimate derived from the same lines. The finding of such a high mass in a face-on object is a strong 
indication that the gas motion is not confined to the accretion disk. The new mass estimate makes 3C 273 accreting with an 
accretion rate about six times lower than the Eddington rate. We discuss the implications of our result for the broad-line-region 
size and black-hole mass vs luminosity relationships for the set of objects for which reverberation black-hole masses have been 
obtained. We find that, while objects with super-Eddington luminosities might theoretically be possible, their existence is not 
necessarily implied by this sample. 
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1. Introduction 

As soon as the variability of the broad emission lines in Seyfert 
1 galaxies and QSOs had been discovered, its potential for 
the investigati on of the s tructur e of these objects became ev- 
ident. Gaskell & Sparke (1986) were the first to propose us- 
ing cross-correlation methods between the continuum and line 
light curves to measure the size of the broad-line region (BLR). 
The subsequent building of long and dense spectro-photometric 
light curves of broad-line active galactic nuclei (BL-AGN), in 
particular tha nks to the IUE sa tellite and to the AGN Watch 
consortium jAlloin et alJ 1 1 9941) . made the reliable measure- 
ment of the BLR size in a significant number of objects pos- 
sible. The so-called reverberation method was later developed 
to allow a dete rmination o f the mass of the cen tral black hole 
(Peterson & Wandel fl999t IWandel et al] fl999). under the as- 
sumption that the gas dynamics is dominated by gravity and us- 
ing the fact that its velocity dispersion can be deduced from the 
width of the line. The connection between black-hole mass and 
line widths in BL-AGN has allowed the development of meth- 
ods to determine the black-hole masses in h igh-redshift QSOs 
(Vestereaard 2002; McLure & Jarvis 2002), and thus opened 
the door to cosmological studies of black holes. These meth- 
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ods require, however, a good understanding of the relationship 
between luminosity and BLR size, which must be investigated 
through the use of the reverberation method. 

iKrolikl J200ll) warned that some assumptions made in the 
context of the reverberation method are not sufficiently val- 
idated. A key assumption is that the gas dynamics is dom- 
inated by the gravitationa l pull f rom the central black hole. 
While Peterson & Wandel (1999) showed that the v oc r ~ 1/2 
relationship expected in such case is satisfied in NGC5548 
for different lines, other competing models can also explain 
this dependence. The best support for the assumptions behind 
the reverberation method is the agreement between the rever- 
beration rnasse^_and_fhose_^taine^j^ dis- 
persion (Geb hardt et all2000l) . lMerritt & FerraresdfeoOll) also 
concluded that the reverberation-mass estimates are not more 
strongly biased than stellar velocity-dispersion ones. The exis- 
tence of hidden assumptions nevertheless makes lKrolikl J2OO 1 ) 
conclude that, even if valid, the method is subject to systematic 
uncertainties as large as a factor 3 or more. 

In a previous paper (Palt ani & Turleil 120031 hereafter 
Paper I), we studied in detail the response of the Lya and Civ 
emission lines in 3C 273. In this paper, we extend our study by 
applying the reverberation method to determine the mass of the 
central black hole in this object. We make the assumption that 
the reverber ation method is vali d and apply the methodology 
presented in Wandel et al. (1999) to derive the black-hole mass 
using the strongest broad ultraviolet lines: Lyo- and Civ. We dif- 
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fer with previous works on the treatment of statistical errors and 
propose here a full error propagation, in order to obtain mean- 
ingful confidence intervals on our results. We al so com pare 
our re sults with those in the optical domain from iKaspi et alJ 
(2000). We finally assess the consequence of our results on two 
important relationships in BL-AGN: the broad-line-region size 
vs luminosity and black-hole mass vs luminosity relationships. 

2. Data 

This paper makes use of data from the 3C 273 database hosted 
by the Geneva Observatory 1 and presented in iTiirler et"al] 
(1999). We use the same data set as in Paper I, namely the IUE 
short wavelength (S WP) spectra processed with the I UE Newly 
Extra cted Spectra (INES) software ( Rodrfgue z-Pascual et alJ 
Il999l) . Details of data preparation can be found in Paper I, 
but we summarize below the main properties of the sample for 
completeness. 

Our light curves each consist of 119 observations cover- 
ing 18 years, including a ten-year period with semi-weekly ob- 
servations during two annual observation periods of 3 months, 
separated by two months. From the original 256 SWP obser- 
vations, over- and underexposed observations were discarded, 
and observations taken within 12 hours were averaged. The 
Lya /U216 and Civ A 1549 fluxes were integrated above a con- 
tinuum defined by a straight line fitted to the points in two 
50 A continuum bands on each side of the lines. We integrate 
the total Lya and Civ line fluxes for projected rest-frame ve- 
locities up to ±20 000 km s -1 and +10 000 km s _1 respectively. 
Uncertainties on the line fluxes take into account the uncer- 
tainty in the underlying continuum and the number of indepen- 
dent bins. The Lya emission line is contaminated by Nv /I1240 
emission. In the following, "Lyc+Nv" will be used to designate 
the full line profile around Lyor, which includes the Nv contam- 
ination, and "Lya-Nv" to designate the line profile around Lyo- 
where the range 3500-10000 km s _1 has been linearly interpo- 
lated to get rid of the Nv contamination. 

Figure ^ shows the three main light curves used in this 
study: the 1250-1300A continuum, the integrated Lya+Nv 
emission line, and the integrated Civ emission line. All light 
curves used in this paper are listed in Table 1 of Paper I, which 
is available electronically. They can also be downloaded from 
the 3C 273 database. 

3. Black-hole mass 

We estimate the mass of the black hole in 3C 273 using the re - 
verberation method introduced by IPeterson &Wandeidl999i) . 
The line-emitting gas is assumed to be in virial equilibrium 
around a mass dominated by the black hole, so that the gas dis- 
tance to the central object and its velocity dispersion are both 
related to the mass through the virial theorem. The velocity 
dispersion can be estimated from the line width, and the gas 
distance from the delay between the continuum and line light 
curves, making possible a determination of the mass. In prin- 
ciple, a single line is sufficient. The consistency of mass esti- 

1 Web site: http://obswww.unige.ch/3c273/ 
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Fig. 1. (a) Ultraviolet continuum light curve averaged in the 
1250-1300 A band, (b) Lya+Nv emission line light curve inte- 
grated for velocities up to 20 000 km s . (c) Civ emission line 
light curve integrated for velocities up to 10 000 km s _I . 



mates from different lines is, however, a stringent test of the 
validity of the method. 

In Paper I, we showed that the Lya and Civ emission lines 
show two distinct components. We interpreted the second com- 
ponent as evidence of gas entrainment in the course of an in- 
teraction with the jet. In principle, the presence of a non-virial 
velocity field makes the above method inapplicable; however, 
we showed that this component mostly affects the blue wings 
of the lines, making the red sides (hereafter, the "red profiles") 
suitable for a mass determination. We nevertheless perform the 
analysis both on the full lines and on the red profiles for com- 
parison. We also perform the analysis separately on Lya+Nv 
and Lya-Nv to evaluate the effect of contamination by Nv. 

In this work, we shall make sure to correctly propagate all 
uncertainties in our data into the determination of the mass. 
As the error distributions of the lag and velocity dispersion are 
not at all Gaussian, they cannot be combined analytically to 
determine the confidence intervals on the black-hole mass. We 
thus construct the full black-hole mass error distribution using 
those of the lag and of the velocity dispersion. 

3.1. Distance of the line-emitting gas 

The distance of the line-emitting gas to the source can be esti- 
mated using the delay between the line and ionizing-continuum 
(or its closest approximation) light curves. Delays can be read- 
ily obtained using a cross-correlation analysis. However, as 
the gas fills a vast range of distances, ther e is no unique way 
to deriv e a delay from a cross-correlation. [Koratkar & GaskelB 
Jl991bl) have argued that the centroid of the cross-correlation 
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Fig. 2. Cross-correlation between the UV continuum and the 
red-profile Ly<?+Nv (top) and Civ (bottom) light curves. In both 
panels, the shaded area represents the parts of the correlation 
above the threshold set at 60% of its maximum (i.e. ft = 0.6). 
Only these areas are used in the estimation. The average lag 
is indicated by a vertical line. Five realizations of the FR/RSS 
method are also shown in each panel for illustration. 



function gives the emissivity-weighted average radius of the 
emission region. 

We follow here the method proposed by IPeterson et alJ 
We calculate the cross-correlations between the con- 
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tinuum (in the 1250-1300A band) and the Lyo-+Nv, Lyor-Nv, 
and Civ light curves. We use t he Interpolat ed Cro ss-Correlation 
Function (ICCF) of Gaske ll & Peterson! \ 19871) w ith the im- 
provements introduced bv lWhite & Petersonld 19941) . The delay 
is estimated by calculating the average lag weighted by the cor- 
relation coefficient for all those parts of the cross-correlation 
which exceed a given threshold ft expressed as a fraction of 
the maximum correlation. This is illustrated in Fig. |2] for the 
actual Lya+Nv and Civ correlations. The lag distribution is es- 
timated using the flux-randomizatio n/rand om-subset-selection 
(FR/RSS) method of lPeterson etafl (fl998): Both line and con- 
tinuum fluxes are first randomized according to their uncertain- 
ties, and 119 dates are drawn at random (with possible repeti- 
tions) from the original dates. The two light curves are then 
cross-correlated, and a new lag estimate is produced. The lag 
error distribution is then constructed from the distribution of 
lag estimates. Examples of FR/RSS realizations are shown in 
Fig El 

Because one wants to take into account as much gas as pos- 
sible, the threshold should be set as low as possible. However, 
if taken too low, non-significant parts of the correlation func- 
tion are included in the lag estimate. Considering the flat-top 
correlation functions we observe here, high thresholds seem in- 
appropriate, as they make the result quite dependent on small 



Table 1. Average rest-frame lags in years of the Lya and Civ 
emission lines. The numbers in brackets give the ±lcr, +2cr, 
and ±3cr confidence intervals. 



Line 



Full profile 



Red profile 



Lyo-+Nv 1.19 
Lyor-Nv 1.21 
Civ 1.89 



1+0.21 1+0.50 1+0.79 

1-0.21 1-0.46 1-0.86 

(+0.20 (+0.49 (+0.79 

1-0.20 1-0.44 1-0.82 

(+0.28 (+0.55 

1-0.29 1-0.53 



J +0.80 
-0.89 



0.97 
0.98 
1.75 



1+0.25 1+0.54 1+0.92 

1-0.24 1-0.63 1-0.85 

(+0.27 (+0.60 (+1.00 

1-0.26 1-0.68 1-0.86 

(+0.39 (+0.69 (+1.08 

1-0.37 1-0.80 1-1.49 




1 2 
Rest-frame Lag [ yr ] 



Fig. 3. Lag error distribution for the red profiles of the Lyor-Nv 
(black), Lya+Nv (light grey), and Civ (dark grey) emission 
lines 



details of the correlation functions. We adopt here a threshold 
of 0.6, but we checked that the choice of the threshold has a 
moderate impact of about 30% on the average delay through- 
out the range 0.3 < ft < 1. The error distributions, how- 
ever, become very large and are no longer unimodal as soon 
as ft > 0.8, which shows that the lags start being dominated by 
non-significant details of the correlation function. The average 
lags and the confidence intervals are given in Tabled All lags 
are expressed in the object's rest frame. In both Lya and Civ 
lines, the lags of the full profiles are systematically larger than 
their corresponding lags in the red profiles, consistent with the 
results of Paper I. The individual differences are not statistically 
significant however. The effect of Nv is completely negligible. 
Error distributions of the lags for the red profiles are shown on 
Fig.E 



3.2. Line velocity dispersion 

The velocity dispersion of the emitting gas can be estimated 
using the shape of the line. In the context of mass determi- 
nation, line velocity dispersion is generally determined from 
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Fig. 4. Continuum-subtracted rms profiles of the Lya+Nv (top) 
and Civ (bottom) emission lines. The grey area marks the re- 
gion of the Lya profile contaminated by Nv emission line. 



the full width at ha lf maximum (FWHM) of the rms profile 
JWandel et al.l 19991) . This profile is then used to avoid contam- 
ination from a possible stable narrow-line component. The re- 
lationship between FWHM and gas velocity dispersion relies, 
however, on assumptions about the line shape that are difficult 
to assess. In addition, the measureme nt of the FWHM is very 
sensitive to details of the profile. Fromerth & Melia (2000) pro- 
pose to calculate the squared velocity dispersion directly from 
the second moment of the line profile: 



<v l >= 3 



2\ 



' j>(v)v 2 dv ( J P(v)vdv 
j>(v)dv ~(fP(v)dv) 



(1) 



where P(v) is the amplitude of the rms profile at rest-frame 
velocity v. The only assumption behind this method is the 
isotropy of the emission, which is the origin of the factor 3 
in the above equation. 

Figure|4]shows the two rms profiles of the Lycf+Nv and Civ 
emission lines. The rms profiles are rather complex, presenting 
a narrow core, broad wings and secondary peaks, and are not 
easily represented by a simple analyt ic form. This is th e kind 
of situation to which the method of From erth & Melial (|2000) 
is particularly adapted; in the following, we adopt Eq. Q to 
calculate the velocity dispersion. 

It is not clear how the error distribution of the velocity dis- 
persion can be calculated. iKaspi et alJ ( 120001) use the FWHM 
dispersion measured in individual spectra as an estimate of the 
uncertainty in the FWHM of the rms profile, but there is no 
obvious reason why it would give the correct result. Here we 
follow a different approach: The uncertainty in any point of the 
rms profile can be estimated using the standard estimator of er- 
ror on the rms: <r rms (A) = rms(A)/ y/2(N - 1), where N is the 



Table 2. Average velocity dispersions in kms~' of the Lye 
and Civ emission lines. The numbers in brackets give the +lcr, 
±2cr, and ±3cr confidence intervals. 



Line 



Full profile 



Red profile 



Lya+Nv 10044 {^g |«g |+n 
Lyor-Nv 9 955 
Civ 



1 119 

341 

+179 (+347 1 + 512 
178 1-366 I- 575 

QfMQ (+352 (+667 ( + 1003 
OOt7 i_ J46 1_ 738 1 |129 



QQ41 !+ 532 (+ 996 !+ 1431 

77 1-539 1-1133 1-1909 

Q7fY7 i+ 240 i+ 461 1+ 674 

y 1-240 1- 491 1- 738 



7 795 



1+442 1+ 828 1+1137 
1-436 I- 947 1-1557 




6000 8000 10000 12000 

Rest— frame Velocity Dispersion [ km s _1 ] 

Fig. 5. Velocity-dispersion error distributions of the red pro- 
files (continuous lines) and full profiles (dashed lines) of the 
Lya-Nv (black), Lycr+Nv (light grey), and Civ (dark grey) rms 
profiles. 



number of spectra. However, the rms in different wavelength 
bins will be strongly correlated with each other, because they 
have been obtained using the same spectra. Thus we scale the 
uncertainties on the rms profile by requiring that the average 
uncertainty of the rms profile in the continuum regions is equal 
to the standard deviation of the rms estimates. We can then es- 
timate the velocity-dispersion distribution by randomizing the 
rms profiles a large number of times according to the uncertain- 
ties in each bin before continuum subtraction (so that the error 
on the continuum determination is included in our determina- 
tion), and by measuring the velocity dispersion using Eq. Q 
for each of them. 

Table |2 gives the velocity dispersions of the Lye and Civ 
emission lines obtained with the 2nd-moment method. Figure|5] 
compares the velocity-dispersion error distribution for the red 
profiles of Lya-Nv and Civ. We removed a spectral bin close to 
the center of the Civ line with a linear interpolation, because it 
was strongly affected by a reseau mark. The effect of removing 
Nv is sufficiently small so that we can exclude a strong residual 
Nv contamination in Lya-Nv. The effect of asymmetry on the 
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Table 3. Average black-hole masses of 3C273 in units of 
10 9 M . The numbers in brackets give the + lcr, ±2<x, and +3<x 
confidence intervals. M is the maximum-likelihood mass ob- 
tained from the combination of Lya-Nv and Civ emission 
lines. 



Line 



Full profile 



Red profile 



Lycr+Nv 
Lya-Nv 
Civ 



.56 
.56 
,.75 



1+1.72 
(-1.68 
i+1.47 
\-1.40 
( + 1.53 
1-1.53 



+4.01 

-3.50 
+3.51 
-3.19 
+3.15 
-2.79 



+6.40 
-6.04 
+5.91 
-5.74 
+4.90 
-4.41 



6.83 
6.59 
7.62 



1 + 1.91 

1-1.88 
(+1.83 
1-1.77 
( + 1.90 
1-1.80 



+4.25 
-4.52 

+4.07 
-4.55 
+3.77 
-3.64 



+7.10 
-6.02 
+6.67 
-5.79 
+6.20 
-6.50 



M 



8.20 



+ 1.02 1+2.31 1+3.82 
-1.14 1-2.25 1-3.16 



6.59 



1 + 1.68 1+3.04 1+4.62 
1-0.90 1-1.87 1-3.15 



measurement of the velocity dispersion is approximately 3% 
in both Lya-Nv and Civ. This consistency shows that the line 
is dominated by a symmetric velocity component, assumed to 
be due to Keplerian motion, but the broadening of the full pro- 
file with respect to the red profile, although not significant, is 
systematic. 



3.3. Mass determination 

The velocity dispersion and delay can be combined to obtain 
the black hole mass from the virial equation: 



M BH = 71.31 t < v > M F 



(2) 



where r is the lag expressed in units of years and < v 2 > the 
squared velocity dispersion expre ssed in units of km s -1 . This 
formula is identical to Eq. (6) of Wandel et al. ( 1999), except 
th at the velocity dispersion is used instead of the FWHM. As 
in Peterson & Wandel ( 1999), the factor between the kinetic 
and potential energies in the virial theorem is set more or less 
arbitrarily to -1. Error distributions of the masses are con- 
structed using a Monte-Carlo approach. We draw a large num- 
ber of pairs of random values from the error distributions of the 
lags and velocity dispersions, and combine them using Eq. Q 
to obtain the error distributions of the masses. The results are 
given in Table[3] As expected from Paper I, the average masses 
obtained with the red profiles of the lines are systematically 
smaller than those obtained with the full profiles. The full- 
profile values fall however within, or close to, the + l<x red- 
profile confidence interval. The difference between Lya+Nv 
and Lya-Nv is less than 5% for the red profiles. The minimum 
mass at the 3<x level is about 0.8 10 9 M . 

The black-hole mass distributions for the red profiles of the 
Lya+Nv, Lya-Nv, and Civ emission lines are shown in Fig. [6] 
All three distributions overlap very significantly. Confronted 
with several estimates of the same values, we can combine 
them using a x 2 test only if the error on the estimates have 
Gaussian distributions. As this is not the case here, we use a 
maximum-likelihood estimator (MLE), where mass is the un- 
known parameter that we try to recover. In the simple case of 
N independent measurements of the same value that we con- 
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Fig. 6. Top panel: Black-hole mass distributions for the red 
profiles of the Lya-Nv (black), Lya+Nv (dashed line), and Civ 
(grey) emission lines. Bottom panel: Likelihood function of the 
mass. The horizontal grey lines show, from top to bottom, the 
±lcr, ±2cr, and +3<x levels. The vertical lines show the M±\, 
M±2, and M ± 3 intervals. 

sider here, the (natural) logarithm of the likelihood function is 
a function of the mass M given by: 



log L(M) = ]T log P t (M), 



(3) 



where P;(M) is the probability distribution of M in measure- 
ment i and is approximated by our black-hole mass distribu- 
tions. The MLE M is then the value of the mass that maximizes 
L{M): 



d log L(M) 
dM 



0. 



(4) 



The {M-„, M+„) n-sigma confidence interval can be determined 
using: 



logL(M ±n ) = logL(M)--. 



(5) 



Applying this method to measurements of the black hole 
mass for the red profiles of the Lya-Nv and Civ emission lines, 
we obtain a maximum-likelihood mass of 6.59 10 9 M . The 
lcr, 2<r, and 3cr confidence intervals are 5.69 - 8.27 10 9 M , 
4.72 - 9.63 10 9 M , and 3.44 - 1 1 .21 10 9 M Q respectively. The 
difference in the MLEs for the full and red profiles is about 

1.5cr. 
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Table 4. Average rest-frame lags in years of the Ha, H/3, and 
Hy emission lines found in this study compared to the val- 
ues of K00. The numbers in brackets give the ±lcr, +2cr, and 
±3cr confidence intervals. Only the ±lcr confidence intervals 
are available from K00. 



Line 


UV-continuum lag 
(This work) 


Optical-continuum lag 
(K00) 


Ha 
Hfi 
Hy 


7 68 1+0.14 1+0.33 (+0.58 
z " uo 1-0.15 1-0.36 1-0.56 
9 ™ (+0.19 (+0.41 (+0.71 
Z ' J0 1-0.19 1-0.33 1-1.20 
9 oc f+0.32 (+0.71 i+1.00 
^■ OJ 1-0.32 1-0.69 1-2.63 


0.90 {-g;|I 
0.77 



4. Mass derived from Balmer lines 

Kas pi et all {2000) (hereafter K00) performed a series of 
spectro-photometric observations in the 4500-8000 A domain 
on a set of 28 QSOs, among which is 3C 273. From the Balmer 
lines, they obtained a mass of 0.235^ ^3 10 9 M from the rms 
profiles, and 0.550^J [^ 10 9 M Q from the average line profiles. 
These values are apparently the average of the masses obtained 
using Ha and H/3. Individual estimates of the mass for different 
lines and different line-width determination methods vary by a 
factor 5 from 0.115 to 0.56 10 9 M , i.e. a much larger range 
than the ~ 15% uncertainties they claim. This indicates a seri- 
ous problem in their mass determination for this object. These 
masses are in complete disagreement with our result using Lya 
and Civ, as they are between 12 and 60 times smaller than 
our maximum-likelihood determination, and at least 6 times 
smaller than M_3. 

It must first be noted that K00 determined the line lag 
with resp ect to con tinuum using the optical continuum around 
5000 A. Palta ni et alJ (ll998) showed that the optical continuum 
properties of 3C 273 differ considerably from those of the ultra- 
violet continuum. They concluded that the optical emission of 
3C 273 is strongly contaminated, and perhaps even dominated, 
by non-thermal radiation, possibly related to the jet of 3C 273, 
and varying on much longer time scales than the UV contin- 
uum. The optical continuum therefore appears unsuitable for 
studying the lag between the ionizing continuum and the lines. 
Using the (full-line) Ha, H/3 and Hy light curves from K00, we 
redetermine the lags with the same ultraviolet continuum light 
curve used to determine the lags with Lya and Civ. In addition 
to using a continuum much closer to the ionizing continuum, 
the lag determination is improved by mutiplying by three the 
number of observations over a twice as long period of time. 
Fig. [7] shows the correlation between the UV continuum and 
the three Balmer lines Ha, H/3, and Hy. Table |4] compares our 
lags with those from K00, ours being between 2 and 4 times 
larger. Furthermore, our estimates are in excellent agreement 
with each other, while in K00 the Ha lag is 60% larger than 
that of Hy. 

The velocity dispersion measurements in K00 are also 
problematic for 3C273. The rms widths of the Balmer lines in 
3C 273 reported by K00 are completely incompatible with each 
other, differing by a factor 2.5, while the quoted uncertainties 
on the widths are only about 2-6%. Such disagreements are ab- 
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Fig. 7. Cross-correlation between the UV continuum and the 
Ha (top), H/3 (center), and Hy (bottom) light curves. In all three 
panels, the shaded area represents the parts of the correlation 
above the threshold set at 60% of its maximum (i.e., § = 0.6). 
Only these areas are used in the weighted average of the lag 
(indicated by a vertical line). The dashed vertical lines indicate 
the lags from K00. Five realizations of the FR/RSS method are 
also shown in each panel for illustration. 



sent in most objects of their sample. As we do not have access 
to the spectra used in K00, we use two spectro-photometric 
observations performed with the STIS instrument (044301010 
for Hy and 044301020 for H/3 and Ha) on the Hubble Space 
Telescope. Without variability information, we have to apply 
the second-moment method to the average profile (calculated 
with a single profile) instead of the rms profile. The contin- 
uum and line regions were defined exactly as in K00. In theory, 
the line width can be reduced by the presence of narrow-line 
Balmer emission, but K00 showed that the effect over their 
whole sample is small. Because of strong [Om] 4959A and 
5007 A contamination in the red profile of H/3, we used only the 
blue profile of this line. Some contaminating emission lines are 
present in the Ha and Hy profiles, but we checked that the blue 
and red profiles give results in good agreement, showing that 
the contamination is probably small. The errors on these values 
were arbitrarily set to 10%, because, if calculated as in the case 
of the UV lines, the formal errors would be unrealistically low. 
Table |3 compares our velocity dispersions with those in K00. 
Except for Ha, which is 80% larger, the widths obtained here 
are comparable to those in K00 using the mean profile. The 
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Table 5. Average velocity dispersions in km s _1 of the three 
Balmer lines found in this study compared to the values of 
K00 with the rms and average (mean) profiles. The numbers in 
brackets give the ±l<x confidence interval. A 10% uncertainty 
has been assumed for our measurements. 
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single profile 


rms profile 
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mean profile 
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Table 6. Average black-hole masses of 3C273 in units of 
10 9 M for the three Balmer lines found in this study compared 
to K00 using the rms and average (mean) profiles. The numbers 
in brackets give the ±l<x, +2<x, and ±3<x confidence intervals. 
The uncertainties on the individual K00 values are of the order 
of 30%. M is the maximum-likelihood masses obtained from 
the combinations of the three Balmer emission lines and of the 
two UV and three Balmer emission lines, respectively. 



Line 
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single profile 


rms profile 


K00 

mean profile 
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HjS 
Hy 
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+2.59 
-1.98 
+ 1.35 
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+2.04 
-2.32 


0.115 
0.363 
0.459 


0.513 
0.564 
0.549 


M (Balmer) 
M (All) 


9 49 /+0-3 4 1+0.66 
1-0.29 1-0.56 

2 44 l+° 51 !+ - 87 
(-0.30 1-0.54 


+1.03 
-0.84 
+ 1.24 
-0.80 







rms velocity-dispersion measurements of K00 for H/3 and Hy 
are about 15-25% smaller, and the discrepancy reaches a fac- 
tor 4 for Ho-. Considering that the K00 values for the different 
Balmer lines using the mean profile differ by 25% and those 
using the rms profile by a factor 2.5, it seems to us that there 
could be a problem in their rms profile determination, which 
could easily result from one or two bad spectra. 

The masses obtained for the three Balmer lines are re- 
ported in Table |6] and Figure [8] shows the mass error distribu- 
tions obtained for the five emission lines. The three Balmer- 
line mass estimates agree rather well with each other, but, 
even though the overlap is not completely negligible, these 
estimates seem significantly lower than the UV ones. The 
lower mass obtained from the Balmer lines can be due to the 
different systematics involved in the UV and optical studies 
and, in particular, to the presence of narrow Balmer emission, 
which artificially decreases the measured velocity dispersion. 
A maximum-likelihood estimate of the mass can be obtained 
by combining the Lya-Nv, Civ, Ha, H/3, and Hy black-hole 
mass error distributions. We obtain a mass of 2.44 10 9 M . The 
lcr, 2cr, and 3cr confidence intervals are 2.14 - 2.95 10 9 M , 
1.90 - 3.31 10 9 M , and 1.64 - 3.68 10 9 M Q respectively. We 
note that the weights of the Balmer lines overwhelm those of 
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Fig. 8. Top panel: Black-hole mass error distributions for the 
five emission lines. The continuous black and dark grey lines 
are, respectively, the Lya-Nv and Civ histograms as in Fig. [6] 
The dashed lines are Ha (black), H/3 (dark grey), and Hy (light 
grey) respectively. Bottom panel: Likelihood function of the 
mass. The horizontal grey lines show, from top to bottom, the 
+ lcr, +2cr, and +3<x levels. The vertical lines show the M ± \, 
M±2, and M ± 3 intervals. 

the UV lines. This is because of their much narrower uncer- 
tainty distributions. However, the UV lines should be less af- 
fected by systematics (the narrow-line emission in particular). 
This maximum-likelihood estimate allows us to set a very se- 
cure 3<x lower limit on the mass of the black hole of 1.6 10 9 
M . Using only the Balmer lines, we find a mass of 2.42 10 9 
M , with lcr,2cr, and 3<x confidence intervals of2. 13-2.76 10 9 
Mq, 1.86 - 3.08 10 9 M , and 1.58 - 3.45 10 9 M , respectively. 



5. Discussion 

5.1. Black-hole mass and Eddington ratio 

The main result of this study is the discrepancy between the 
masses found by K00 and our estimates, since they differ by a 
factor ranging from 12 up to 60. As our study of the Balmer 
lines shows, discrepancies affect both the estimates of the de- 
lays and of the line widths. A general argument in favor of 
our mass estimate is the good qualitative agreement between 
the individual estimates for the Lya, Civ, and the three Balmer 
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lines. The largest difference between the UV and Balmer lines 
is a factor three, while the systematics in the two sets of data 
can be quite different. By contrast, the K00 estimates, reported 
in Table |6] differ by a factor 5 between each other, and about 
4.5 for Ha only, although all the estimates have been obtained 
with a consistent data set. In our study, we could not deter- 
mine the rms profiles of the Balmer lines (not even their aver- 
age profiles), because of the lack of data. Therefore, it is quite 
probable that a large part of the UV-Balmer discrepancy would 
be alleviated if the same methodology could be applied to the 
Balmer lines as well. It must be noted that all Balmer profiles 
are affected by contaminating emission lines, in addition to the 
narrow-line Balmer emission: [Nn] /16548 and [Nn] /16583 for 
Ha; [Ora] ,14959 and [Om] A5007 for H/3; [Om] ,14363 for 
Hy. The [Om] lines being particularly strong in the wing of 
H/3, its red profile had to be excluded from our analysis. The 
HjS region is moreov er seriously c ontaminated by Fen pseudo- 
continuum emission (B oroson & Greenll992l) . making the line 
profile of H/3 quite uncertain. The remaining narrow lines, plus 
the narrow-line Balmer emission, may significantly modify the 
true broad-line profiles. As they fall close to the Balmer line 
centers, their effect is to reduce the line second moments, pro- 
viding underestimated black-hole masses. We thus conclude 
that, although affected by larger statistical errors, the UV-line 
estimates very probably suffer from smaller systematic errors. 
Therefore, the best mass estimate is the combined UV mass: 
6.59 10 9 M , with a 3cr lower limit of 3.44 10 9 M Q . Including 
the Balmer lines, the 3<x lower limit is found to be 1.6 10 9 M . 
This lower limit is very robus, and should only increase if the 
systematic errors are reduced. 

We point out that this fairly high mass estimate is ob- 
tained under the hypothesis that the gas motion is virialized 
and (close to) isotr opic. As 3C273 has a superluminal jet (e.g., 
ICourvoisierl fl998), 3C273's accretion disk is very probably 
close to face-on (assuming the disk is perpendicular to the jet), 
so projected velocities of particles in the disk would be greatly 
reduced by the geometry. If the emission lines originate from 
the disk, the mass estimate should be multiplied by a factor 
1 / sin i > 6, if the angle between the perpendicular to the disk 
and the line of sight is i < 10°. On the basis of the response 
of the UV lines, we also concluded in Paper I that accretion- 
disk geometries of the broad line region can be excluded in this 
object. 

Our comparison between the red and full profiles for the 
Lyo- and Civ emission lines shows that inclusion of the possibly 
outflowing gas seen in Paper I would misleadingly increase the 
estimate of the mass of the black hole by approximately 15- 
30%. This is far from negligible, but it also shows that the order 
of magnitude of our result is not affected by the choice of using 
only the red profile. 

From Courvoisier (1998) and lTiirler et al.l \ 19991) . we esti- 
mate the bolometric luminosity of 3C 273 to be approximately 
lO^ergs- 1 (using H = 70^8-' Mpc" 1 ), giving an Eddington 
mass of about 1O 9 M . Using the 5 100 A flux as an estimate 
of the bolometric luminosity and a different cosmology, K00 
find a comparable value of ~ 0.5 10 9 M o . According to the 
results of K00, the ratio between the mass derived from the 
Eddington luminosity and that from the reverberation method 



is therefore in the range l-A, depending on which mass esti- 
mates are used, which makes 3C273 a very strong candidate 
for super-Eddington accretion. Using our new reverberation- 
method mass, the Eddington ratio is estimated to be about 0.15 
and, using the individual estimates, falls in the range 0. 13-0.53, 
i.e. in a fully sub-Eddington regime. 

5.2. Size of the broad-line region 

The size of the Balmer-emitting broad-line region in AGN 
has been investigated by s everal authors (Koratkar & Gaskell 
Il991al IWandeletalJll999l K00). While the first two studies 
find a BLR size Rblr scaling roughly as L , K00 find a steeper 
relationship Rblr ~ L°- 700±0 033 using a sample of 34 objects 
that were investigated using the reverberation method. The L 5 
relationship is expected if gas density and ionization parame- 
ters are independ ent of the source lu minosity in the emitting 
region. However, Vestergaard (2002) correctly points out that 
the K00 relationship is calculated with a method unable to 
take the strong intrinsic scatter into account, which produces 
biased slo pes and too low uncertaint ies. Applying the BCES 
method of lAkritas & Bershadvl (fl996) on the K 00 data and us- 
ing a modified BLR size for NGC4051 from IPeterson et all 
(2000]), she found Rblr ~ L () 58±0 09 , perfectly compatible with 
L . We can refine this relationship by using the new Rblr 
for 3C273 found in t his work and by includin g new measure- 
ments for NGC 378 3 JOnken & Peter sonl2002l) . and NGC 3227 
JOnkenetalJ2 003). 3C 273 has an important weight in the rela- 
tionship, because it is the highest-luminosity object in the sam- 
ple. From a maximum-likelihood estimator on the lag of the 
three Balmer lines, we have Rblr = 986*2J days. Applying the 
BCES estimator, we obtain 

/?blr~^(5100A)°- 59±0 09 . (6) 

The two relationships are perfectly compatible with each 
other and still compatible with a L 5 relations hip. Excluding 
NGC 4051 as proposed in Vestergaard (2002), the slope be- 
comes 0.67 + 0.10, but it is not clear why this object should 
be discarded, as its distance to the BCES estimator does not ap- 
pear unreasonable, considering its uncertainty. It must be noted, 
however, that this object has a very strong weight in the rela- 
tionships 

iMcLure & Jarvisl J2002I) recalculated the ^(51 00 A) lumi- 
nosities using fluxes from lNeugebauer et al 1 11 19871) and a mod- 
ern cosmology (Q. m = 0.3, Qa = 0.7, Hq = 70kms~' Mpc -1 ), 
and found a slope flatter than K00. We find here that, while 
the slope of the bisector is slightly lower, it is not significantly 
affected by the new luminosities, showing that the K00 lumi- 
nosities are sufficiently accurate. 

5.3. Mass-luminosity relationship 

With the same sample, K00 investigated the mass-luminosity 
relationship in 34 Seyfert 1 and QSOs. They found that the re- 
lationship deviates significantly from the L ~ M relationship, 
which is expected if the Eddington ratio Lboi/^Eddington is in- 
dependent of the black-hole mass. They find Mbh ~ L°- 5±01 , 
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Fig. 9. Mass-luminosity for the 34 objects in the K00 sam- 
ple. The solid line is the BCES bisector. The dashed line is 
the BCES bisector obtained after discarding unreliable objects. 
The shaded area shows the Icr error domain. The square, circle, 
triangle and pentagon respectively show 3C273, NGC 4051, 
PG 1613+658, and PG 1704+608. 



meaning that Seyfert 1 galaxies have low accretion rates com- 
pared to QSOs. It even strongly suggests that very luminous 
QSOs must accrete in the super-Eddington regime. This result 
has indeed sparked a lot of intere st in the theoretica l investiga- 
tion of super-Eddington models dCollin et all2002l) . 

As we find that 3C 273 has a lower accretion rate than pre- 
viously thought, we reexplore the mass-luminosity relationship 
using the BCES bisector. With the new 3C 273 mass found here 
and the updated masses of the three Seyfert galaxies, we find 



M BH ~ AL A (5 100 A) 



0.65±0.18 



(7) 



Figure [5] shows the updated mass-luminosity relationship. 

We caution, however, that this relationship can be disturbed 
by a few ill-measured masses. For instance, in the K00 sample 
two objects, PG 1613+658 and PG 1704+608, have masses dif- 
fering by a factor 5 or larger, depending on whether the rms- or 
mean-profile method is used. If we discard these objects as too 
unreliable, we obtain 



6. Conclusion 

The most important consequence of this work is to show that 
the reverberation method can provide quite different results 
on a single object, as we find a black hole mass of about 
6.59 10 9 M , an order of magnitude larger than in the previous 
study o f K00 . This result is considerably stronger than that of 
Krolikl J200lh . because it shows that serious discrepancies can 
result from the measurement of the delays and line widths and 
not only from the underlying assumptions of the method. In this 
study, we benefited from a much better determination of the de- 
lay between the lines and the continuum, thanks to the use of 
the ultraviolet continuum with a better sampling. However a 
large fraction of the discrepancy is caused by the determina- 
tion of the line widths, for which we can propose only tentative 
explanations. 

The mass determined from Balmer lines is about a factor 
3 smaller than that obtained from the UV lines, but because 
of the lack of adequate data we could not determine the rms 
profile of the Balmer lines and had to rely on a less optimal 
velocity-dispersion measurement method. Thus at least part of 
the discrepancy can be explained by the different systematics 
and in particular by the inclusion of narrow-line emission in 
the line profile. Our Balmer-line mass is, in any case, at least a 
factor 5 larger than the previous study's estimate. We are able 
to put a very robust 3<x lower limit on the mass of the black 
hole in 3C273 of 1.61 9 M . 

We studied the mass-luminosity relationship of supermas- 
sive black holes determined with the reverberation method. 
Using our new 3C273 mass estimat, and updated estimates for 
three Seyfert galaxies, we find that the case for a correlation 
between the Eddington ratio and the black-hole mass in the 
reverberation-method sample is far from established. 

We add finally that, in a very recent in-depth reanalysis 
of the exhaustive data set from the AGN Watch collaboration 
jAlloin et alJll994h . lPeterson et ail J2004 obtained a new esti- 
mate of the mass of 3C273 of 0.89+0.19 10 9 M o . Using our de- 
lay estimates for the Balmer lines, the estimated mass would be 
approximately 2.5 1O 9 M , perfectly compatible with our own 
estima te using the Balmer lines (2.42 10 9 Mq). IPeterson et alJ 
(2004) also give a much improved estimate of the slope of the 
mass-luminosity relationship Mbh ~ L°- 79±0 09 , again in com- 
plete agreement with our work. 
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M BH ~ ALa(5 100 A) 



0.79±0.19 



(8) 



The reverberation-mass sample is therefore compatible 
with an Eddington ratio independent of the black-hole mass . 
This is consistent with the conclusion of I Woo & Urrvl d2002t) . 
who investigated the mass-luminosity relationship in a sample 
of 377 AGN whose masses had been estimated using different 
methods, including reverberation. If the absence of relationship 
is confirmed, it means that there may currently be no need for 
super-Eddington-luminosity models. 

He re as well, using the luminosities from lMcLure & Jarvisl 
(2002i) does not significantly affect the bisector estimates. 
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